Program LMTART for Electronic Structure Calculations 
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A computer program LMTART for electronic structure calculations using full potential linear 
mufBn-tin orbital method is described. 
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I. DESCRIPTION 

Full-potential linear-muffin-tin-orbital (Ref.^) (FP- 
LMTO) program LMTART is designed to perform band 
structure, total energy and force calculations of solids us- 
ing methods of density functional theory (DFT) (Refi^). 
The development of LMTART has been initiated in 
1986 in P.N. Lebedev Physical Institute. It was fur- 
ther extensively developed in Max-Planck Institut fiir 
Festkorperforschung, and its most recent contributions 
have been done in Departments of Physics of Rutgers 
University and of New Jersey Institute of Technology. 

LMTART performs electronic structure calculations 
using popular local density approximation (LDA) and 
generalized gradient approximation (GGA) of density 
functional theory^. To deal with strongly correlated sys- 
tems it also implements LDA-I-U method- and will in- 
clude the capability to study materials using spectral 
density functional approach and dynamical mean field 
theory (DMFT)4 in the nearest future. For determining 
equilibrium crystal structure and phonon spectra LM- 
TART computes total energies, forces^ as well as full 
wavevector dependent lattice dynamical using linear re- 
sponse theory. For studying optical spectral functions 
LMTART evaluates dipole matrix elements between var- 
ious electronic states and can determine optical conduc- 
tivity of materials. 

The code implements atom centered local mufSn-tin 
orbitals as a basis for representing one-electron wave 
functions in both bare and screened (tight-binding) rep- 
resentations with general choice of tail energies k^. This 
allows to increase accuracy in the calculations by using 
multiple-K LMTO expansions. Low lying core electronic 
states can be resolved as bands in separate energy panels. 
For heavy elements, spin-orbit coupling matrix elements 
can be taken into account using a variational procedure 
proposed by Anderseni. Full three dimensional treat- 
ment of magnetization in relativistic calculations includ- 
ing LDA-f U is implemented to study problems related to 
magnetic anisotropy and non-coUinear magnetism. Ef- 
fect of finite temperature can be modelled by introducing 
Fermi statistics for the one-electron states. Additionally, 
LMTART can determine tight-binding fits to the energy 
bands by extracting hopping integrals between various 
orbitals using tight-binding LMTO representation. 

LMTART is written on FORTRAN 90 and uses dy- 
namical memory scheme. No additional recompilation of 
the code is required when changing numbers of atoms, 



spins, plane waves, etc. There are two basic regimes for 
working with this program: the self-consistent charge 
density calculation and calculation of physical properties 
such as electronic structure, optical properties, etc.. The 
simplest input to LMTART involves only atomic charges 
of the atoms as well as crystal structure. 

LMTART works with two different approximations re- 
lated to a shape of the potential: (i) atomic sphere ap- 
proximation (ASA) and plane wave expansion (PLW). 
ASA uses overlapping atomic spheres, where the poten- 
tial is expanded in spherical harmonics inside the spheres, 
but any contribution from the interstitial region is ne- 
glected. Such method is fast and provides reasonably 
good energy bands. However it is not sufficiently accu- 
rate to deal with distortions and phonons. PLW is a 
full potential approximation which uses non-overlapping 
muffin-tin spheres, where the potential is represented 
via spherical harmonics expansions, and the interstitial 
region where the potential is expanded in plane waves. 
The full potential regime provides the best accuracy at 
the price of increasing computational time. Finally, LM- 
TART can be run in a special tight-binding regime by 
supplying hopping integrals between various orbitals. 

Two extensions of LMTART are available as separate 
programs. One is a full-potential linear-response linear- 
mufhn-tin-orbital package LMTO PHONONS which is 
designed to perform linear-response calculations of the 
phonon spectra for arbitrary wave vectors q. Main fea- 
tures of this code include: (i) Computations of first-order 
changes in the charge density and the potential due to 
displacements of nuclei for arbitrary wave vectors q. (ii) 
Calculations of the dynamical matrix and the phonon 
spectra, (iii) Calculations of elect ron-phonon interac- 
tions, Eliashberg spectral functions a^Fiuj) and trans- 
port properties such as electrical and thermal resistivi- 
ties. The description of the method and complete set of 
references can be found in Ref.^. Numerous applications 
are given in Ref.^. 

Also available is an extension of LMTART to compute 
dynamical susceptibilities and magnon spectra This is 
a full-potential linear-response linear-muffin-tin-orbital 
package of programs LMTO MAGNONS designed for 
these purposes. As a by-product, spin wave spectra 
and their lifetimes are accessible as peaks in calculated 
imaginary spin susceptibility. A short description of the 
method can be found in Refi^. 

Software MINDLab which runs under Microsoft Win- 
dows operating systems can be used to set up input files 
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and analyze output files of LMTART. It can also be 
used to run the code without learning its extensive in- 
put options. This software is necessary to visualize all 
data withdrawn from LMTART. Currently, MINDLab 
can perforin the following tasks (i) Crystal group cal- 
culations. By setting atomic positions MINDLab shows 
crystal group operations found for given atomic config- 
uration; (ii) Fat bands calculations. Fully colored visu- 
alization of separated orbital characters on top of the 
band structure can be performed by mouse click oper- 
ations, (iii) Density of states calculations. MINDLab 
computes and visualizes densities of states, total and or- 
bital resolved, (iv) Extractions of hopping integrals for 
tight -binding fits can be performed quickly with MIND- 
Lab. (v) Optical properties calculations. Dielectric func- 
tions ei(ti^), £2{^), and electron energy loss spectra can 
be computed and visualized by MINDLab. (vi) Visu- 
alizations of charge densities, full potentials, and Fermi 
surfaces in both two and three dimensions using OpenGL 
graphical library implemented within MINDLab (viii) 
Three dimensional visualization of crystal structures and 
their distortions, (ix) Full three dimensional visual- 
ization of vector fields including velocity fields on top 
of the Fermi surfaces and magnetizations, (x) Simplified 
input for studying correlated electronic systems which 
currently implements LDA-f U method and will include 
LDA-f DMFT technique in the nearest future. 

MINDLab works with LMTART as with an executable 
file which reads input data, performs band structure cal- 
culation, and stores output files. MINDLab controls 
this process, it prepares input for LMTART using dialog 
windows, and executes LMTART program as a separate 
thread. When the calculation is finished, LMTART noti- 
fies the MINDLab, and the output files can be visualized 
by simple mouse click operations. While running LM- 
TART for complicated compounds can be rather slow 
(it depends how many atoms per unit cell is selected), 
MINDLab can be used to prepare the input, after which 
the files can be copied to a more powerful computer where 
the LMTART can be executed much faster. After run is 
performed, the output can be copied back to the PC and 
analyzed quickly by MINDLab. All input/output files of 
LMTART are formatted ASCII files, and therefore sys- 
tem independent. 



II. APPLICATIONS 

Several most recent applications of LMTART can be 
found in Refs4*2*2i2iiS. They include calculations of 
phonons and electron-phonon interactions in novel su- 
perconductors, such as MgB2,LiBC, MgCNis, studies of 
phonons in strongly correlated systems, e.g. NiO, Pu, 
computations of magnetic anisotropy energies of ferro- 
magnets Fe, Co, Ni, Cr02 evaluations of dynamical spin 
susceptibilities and spin wave spectra, and so on. 

Here we illustrate one of the most recent calculationii 
of phonon spectrum in MgCNis using a linear response 



regime of LMTART. The calculated phonon dispersions 
along major high symmetry lines of the cubic Brillouin 
zone are given on Fig. 1. The frequencies are seen to 
be span up to 1000 K, with some of the modes showing 
significant dispersion. In general, we distinguished three 
panels where the top three branches around 900 K are 
carbon based, the middle three branches around 600 K 
are Mg based and 9 lower branches are all Ni based. 
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FIG. 1: Calculated phonon spectrum of MgCNis along major 
symmetry directions of the Brillouin zone of the cubic lattice 
using density functional linear response method. Some curves 
are widened proportionally to the phonon linewidths. 

We have discovered a striking feature of this phonon 
spectrum connected to the presence of a low-frequency 
acoustic mode which is very soft and is even seen to be 
unstable along some points in (^^0) direction in the Bril- 
louin zone. This mode is essentially Ni based and cor- 
responds to perpendicular movements of two Ni atoms 
towards octahedral interstitials of the perovskite struc- 
ture. The latter is made of each of the four Ni atoms and 
two Mg atoms. 

The persistence of the instability which does not oc- 
cur for any of the high-symmetry point needs a non- 
trivial frozen-phonon analysis. As our polarization vec- 
tors prompt that the largest anharmonicity is expected 
for the rXM plane (qz = 0) of the Brillouin zone of the 
cubic lattice, we have performed three such calculations 
for the points a =(ii0)^, b =(iiO)^, and M. The re- 
sults of these calculations reveal essentially anharmonic 
interatomic Ni potentials. A shallow double well with a 
depth of the order of 20-40 K and the curvature at the 
equilibrium of the order of 40~50i K exists for the a and b 
points which becomes vanishing at the M point. We have 
found such a behavior by both total energy and force cal- 
culations for the supercells of 20 atoms (a and b points) 
and 10 atoms {M point). Clearly, such a small depth 
on the temperature scale indicates that the distortions 
are dynamical and zero point motions would prevent the 
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appearance of the static long-range order. Complete dis- 
cussion of these calculations as well as the comparisons 
with existing experiments can be found in Ref.^. 

Fig.l also shows the calculated phonon lincwidths 
7qj,by widening some representative phonon dispersion 
curves ujqi, proportionally to j^^. In particular, we 
found that some phonons have rather large linewidths. 
This, for example holds, for all carbon based higher ly- 
ing vibrational modes. The strength of the coupling, 
Aq^, for each mode can be obtained by dividing by 
7riV(eF)Wq^, where N{eF) is the density of states at the 
Fermi level equal to 5.3 st./[eV*cell] in our calculation. 
Due to large OJ^^, this unfortunately results in strongly 
suppressed coupling for all carbon modes which would 
favor high critical temperatures. The coupling, however, 
is relatively strong for the Ni based modes. For example, 
we can find A's of the order of 1-3 for the Ni based optical 
phonons around 250 K. Our resulting value of electron- 
phonon coupling constant A is found 1.51. This assumes 
that MgCNis is a strongly coupled electron-phonon su- 
perconductor. 

III. DOWNLOADS 

All LMTART based programs described above are 
available for downloading at the following URL: 



http://www.physics.njit.edu/~mindlah. Each program 
consists of the source code, manual, and example files. 
An extensive database of electronic structures for differ- 
ent materials is also available at this HTTP URL. 
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